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Abstract. Quenched thermodynamic states of an amorphous ferro- 
magnet are studied. The magnet is a countable collection of point par- 
ticles chaotically distributed over R'', d > 2, which is modeled by a 
homogeneous Poisson point process. Each particle bears a real-valued 
spin with symmetric a priori distribution; the spin-spin interaction is 
pairwise and attractive. For every pair of particles, the interaction in- 
tensity is random with distribution dependent on the Euclidean distance 
between the particles. The intensities are independent of the underlying 
Poisson point process and also of each other for distinct pairs of parti- 
cles. For this model, we prove that with probability one: (a) quenched 
thermodynamic states exist; (b) they are multiple if the particle density 
and the interaction strength are large enough. 



1. Introduction 

In this paper, we study thermodynamic states of the following model of 
an interacting particle system. A countable collection of point 'particles' 
is chaotically distributed over M°', d > 2. The corresponding mathematical 
model is a homogeneous Poisson point process vr^ with intensity A > 0. Each 
'particle' represents a cluster of magnetically active physical particles, and 
hence is supposed to bear spin ax which takes any real value. We assume that 
ax is characterized by a symmetric a priori distribution x on M, the same for 
all X. The spin-spin interaction is supposed to be pair-wise and attractive. 
For the 'particles' located at x and y, it has the form JxyfTx^^y with intensity 
Jxy = (pi\^ ~ y\)Dxy, where a non-random (measurable) function (p takes 
values in [(/>=„,(/)*], 0=,, > 0, (p* < oo. The random variables {Dxy : x,y G 
W^, X y} take values 1 and with probability g{\x — y\) and 1 — g{\x — y\), 
respectively. They are mutually independent, and also independent of the 
underlying Poisson point process. We suppose that g{r) G [0, 1], and g{r) = 
whenever r > r^,, where r.,, > is a fixed parameter of the model. The 
physical meaning of the factors Dxy is to take into account that some of the 
exchange interactions between the spins can be suppressed to zero by the 
random environment in which the system is placed. We call this model the 
amorphous ferromagnet, cf. [231 Section 11]. 

In view of the randomness mentioned above, the notion of the thermody- 
namic state of our model can be introduced in the following two ways. In the 
first one, the randomness is taken into account already at the level of local 
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states defined on the space of joint configurations of particles, spins, and 
connection variables D. The global Gibbs states constructed in this way are 
then annealed states; they describe the equilibrium of the whole system. In 
the case of non-random D^y = 1, |x — y| < r^, and D^y = 0, ja; — y| > r^,, the 
mentioned configuration space would be the space of marked particle con- 
figurations 7 = {(x, ax) ■ X G 7}, where 7 is a locally finite subset of M'^, see 
()2.ip below. The second approach, which we follow in this paper, consists in 
constructing thermodynamic states of the spin system alone for fixed typi- 
cal configurations of the particles and the variables D. These are quenched 
states. The global observables characterizing such states are self- averaging, 
i.e., non-random. Note that studying quenched states is a more difficult 
problem, as compared to that of annealed ones, in view of the present spa- 
tial irregularities which do not allow for applying here most of the methods 
effective for regular systems. 

Actually, there are only few publications on the mathematically rigor- 
ous theory of phase transitions in spin systems of general type living on 
non-crystalline (amorphous) substances, cf. [lOl [TTl [121 ES] where annealed 
states were considered. The reason for this is presumably that the meth- 
ods for studying such phenomena, e.g., infrared estimates, are essentially 
based on the translation invariance (and other symmetries) of the under- 
lying crystals. At the same time, for Ising spins ax = ±1, there exist 
methods applicable to the corresponding models on graphs, cf. [THl IT^ [2U]. 
For such models, the main idea of proving the existence of phase transi- 
tions is to relate the appearance of multiple phases of the spin system to the 
Bernoulli bond percolation on the underlying graph. Amorphous substances 
can be described in terms of random point processes, and thereby can also 
be considered as random graphs, in which one can observe a Bernoulli bond 
percolation, see [71[2T1[22]. In the present work, we combine these methods 
and prove that the mean magnetization in the model of an amorphous ferro- 
magnet with spins cia; G M outlined above can be positive almost surely, and 
hence the quenched Gibbs states can be multiple, if the particle density and 
the interaction strength are large enough. We do this as follows. First we 
establish the existence of the corresponding quenched Gibbs states. For ran- 
dom graphs with unbounded vertex degrees, proving the existence of Gibbs 
states with properties suitable for physical applications is a nontrivial prob- 
lem, especially if the single-spin distribution x has noncompact support. In 
the latter case, there can exist states supported on configurations of spins 
with rapidly increasing \ax\ as |x| — t- +c«, whereas for typical ferromagnetic 
configurations in physical substances, most of the spins take values close 
to same s > 0. Therefore, for our model Gibbs measures of physical rele- 
vance ought to be supported on the configurations with tempered growth of 
\ax\- In this paper, the existence of such tempered Gibbs states is proven 
by means of a result of \T6] and a property of the Poisson point process 
obtained in Proposition 12.21 below. Next, by means of results of [71 12H [22]. 
we conclude that the underlying graph based on the Poisson point process 
with the adjacency relation established by the function g almost surely has 
an infinite connected component, in which the Bernoulli bond percolation 
takes place if the intensity A exceeds some threshold A* G (0, 1). By means 
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of a result of , from this we deduce that the Ising model on such a graph 
can be in a ferromagnetic state. Then the generalization to other types 
of the single-spin measures Xi including those corresponding to unbounded 
spins, is performed by means of the Wells inequality [27J. For the reader 
convenience, we present here a complete proof of the latter. 

Finally, let us mention that, for our model with cTx S the problem 
of uniqueness of Gibbs states - inverse to the one which we study here, 
remains open except for some special cases, see subsection 12.41 below. For 
the Ising model, the mentioned uniqueness can be established by showing 
that the Bernoulli site percolation on the underlying graph is absent for 
small enough values of the corresponding probability parameter, see [9]. 

2. Quenched Gibbs measures 

2.1. The model. By P we denote the space of all locally finite subsets 
(configurations) of W^, that is, 

(2.1) r = {7 C M'^ : |7 n ir| < oo for any compact K C M"* }, 

where \A\ stands for the cardinality of a finite set A. This space is equipped 
with the vague topology being the weakest one in which the maps F 9 7 1— > 
'Ylx^'yf^'^) continuous for all continuous functions / : M'^ — )• M with 
compact support, see, e.g., @j for more detail. This allows for introducing 
the corresponding Borel cr-field BiV). The vague topology is metrizable in 
such a way that the corresponding metric space F is complete and separable. 
For A > 0, by tta we denote the homogeneous Poisson measure on (P,;B(P)) 
with intensity A. It is convenient for us to consider ttx as the probability 
distribution of a point process on a complete probability space (r2, J^, P). In 
that we assume the existence of a measurable map O 9 a; 1— )• 7(w) E P such 
that, for each A G B{T), irxiA) = F{-f'^{A)). 

For r=K > 0, let g : [0, +00) — )• [0, 1] be a non-increasing function with 
support in the interval [0, r,,,]. Following [22], we define a system of random 
variables D := {D^y = Dy^, x,y G W^, x y} on {Q,J^,¥) such that each 
Dxy takes values 1 and with probability g{\x — y\) and 1 — g{\x — y\), 
respectively. The function g is assumed to be such that 

(2.2) / gi\x\)dx=:g,£{0 + oo). 

All random variables Dxy are mutually independent and are also indepen- 
dent of the underlying Poisson point process. Now, for a fixed pair (7,!)), 
we consider the graph G(7, D) with vertex set 7 and edge set 

E{-,,D) = {{x,y} Cj:Dxy = l}. 

It is a random graph called the random connection model with connection 
function g, driven by ttx, see [71|22]), and especially [211 pages 18-20] for a 
more detailed exposition of this model. It will serve us as the underlying set 
for spin configurations of the magnet we consider. 

Let B{M.'^) and i3b(K'^) stand for the set of all Borel and all bounded Borel 
subsets of M*^, respectively. For A G B{M.'^) and 7 G P, by Ea(7, D) we denote 
the set of edges of G(7, D) both endpoints of which are in 7a := 7 H A. To 
each X G 7, we assign a spin variable ax G K. Then the configuration of spins 
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corresponding to ^ is cr — {(^x^xf^-y £ M'^. The set of all spin configurations 
My is equipped with the product topology and with the corresponding Borel 
cr-field B{Wf). For A G ^(M"^), we write = M'^ \ A. For such A, by cta 
we denote the 'configuration in A', i.e., cja := : x S 7a}- Given two 
configurations a and a, by cja x we mean the configuration such that 
its restriction to x G 7a (respectively, to x £ ^\c) is ax (respectively, ax)- 
For X S 7, by dx we denote the neighborhood of x in G(7,Z?), that is, 
9x := {y G 7 : Dxy = 1}. 

Let X be a finite symmetric measure on M such that, for some u> 2 and 
X > 0, the following holds 

(2.3) / exp(x|t|")x(dt) < oo. 

Our aim is to construct Gibbs measures on W that correspond to the single- 
spin measures Xx = X ^"^^ to the following relative energy functionals 

(2.4) - E']^{a !s\a f^c) = ^ Jxyf^xf^y + ^ JxyCTx^y 

{x',y}eEA{7,s) a;e7A yefen7Ac 

Here the interaction intensities take the form 

(2.5) Jxy = - y\)Dxy, 

and hence are random. In (|2.5p . <j) : — ?■ M+ is a non-random measurable 
function such that, for some 0* < +oo, 

(2.6) </.(r) G for all r € [0,r,]. 
For A G Sb(IK'^) and a GW,we define 

(2.IJX(^I^) = I Ia{(TA X CTA=)exp(-S2(^A|^A0)XA(dCTA), 

^Alo'j Jrta 

where lyi is the indicator of ^ G ;B(M^), E'^ is as in ()2.4p . and 

(2.8) XA(d^A) := (g) x(dcT,), 

xe7A 

^a(o-) := / exp(-£'2(crA|a-A0)XA(c?crA)- 

JRT 

Thus, for each ^ G B{W), n^(^|-) is i3(]R''')-measurable, and, for each 
a G W^, n^(-|cj) is a probability measure on (M'^, i3(M''')). The collection of 
probability kernels {H^ : A G ;Sb(]R'^)} is called the Gibbs specification of the 
model we consider, see [U Chapter 2]. It enjoys the consistency property 

/ nl^{A\a)Ul^{da\a)=Ul^{A\a), 

which holds for all A G B{W), a G W, and ah Ai,A2 G Bi,{R'^) such that 
Ai C A2. 

Definition 2.1. A probability measure fj, on (M''', i3(M'^)) is said to be a 
quenched Gibbs measure of the model considered if it satisfies the Dobrushin- 
Lanford-Ruelle equation 

^l{A) = [ Ul{A\a)n{da), for all A G B{M.^'). 

JR-t 
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The set of all such measures is denoted by G{'y,D). 

A priori it is not obvious whether Q{'j,D) is nonempty. If this is the 
case, then Q('j,D) depends on the random parameter co through 7 and D, 
not necessarily in a measurable way, cf. [T7j. In this work, we leave aside 
measurability issues, which we are going to address in an accompaining work 

As noted above, the set G{j, D) may contain elements which are not suit- 
able for describing thermodynamic states of a ferromagnet since we have no 
a priori information concerning the support of the eventual ^ G ^(7, D). At 
the same time, Gibbs measures of physical relevance ought to be supported 
on the so called tempered spin configurations, for which \ax\ increases 'not 
too fast' as — -|-cxd, cf. ()2.13p below. Of course, the existence and prop- 
erties of the Gibbs measures depend on the properties of the graph G(7, D), 
which we study in the next subsection. 



2.2. The underlying graph. 



2.2.1. Estimating the degree growth. For x G 7, let n{x) be the number 
of neighbors of x in G(7,Z)), i.e., n[x) := \dx\. Clearly, n{x) is almost 
surely finite since each 7 is almost surely locally finite. Note, however, that 
sup^g^n(x) = +00, also almost surely. 

For an a > 0, we introduce the weight function 

(2.9) Wa{x) = exp(-a|j;|), x G M"'. 

For a; G 7 and ^ > 0, we then consider, cf. Eqs. (4) and (5) in |16j . 

(2.10) a{a,e):= ^ K(x) + u;„(y)] [n(x)n(^/)]^ 

{x,y}€E{j,D) 

h{a) := y^Wgjx). 

Proposition 2.2. For each positive a and 6, both a{a,9) and b(a) are 
almost surely finite. 



Proof. By the very definition of the Poisson measure vr;^, for each n G N and 
any measurable and symmetric function / : M*^ x F — )• M_|_ := [0, +cx)), we 
have that 



(2.11) /(^E/(^'7\x)^ 



TTxid'y) 

^ I { f{x,l)dx ] TTxid'y) 



(the Mecke identity). Then 

Mb{a) = / y^Wa{x) TTxidj) = A / Wa{x)dx < 00, 
Jr J Jwi 
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where E is the expectation with respect to P. Hence b{a) < oo almost surely. 
Next, we rewrite (j2.10p in the form 

(2.12) a{a,0) = 'Ywaix)m0{x), mg{x) := [n{x)n{y)]^ . 

Let X : {0, 1} be the indicator of the bah B2r, := {x E M'^ : |x| < 2r*}. 

Clearly, 

max n{y) < > I{z — x). 
ye{x}Udx ^-^ 

Applying this in (j2.12p we get 

26+1 



meix) < [n(x)]^~'""^ max[n(y)]^ < \^ I{y — x) 
yedx \ ^ 

By this and ()2.1ip . as well as by the translation invariance of vr^, we then 
obtain from (f2J2]l 



Wa{x) ^ I{y - x) 




Here 



261+1 

■K\{d-f) 

26+1 

I 'Kxid-y) ) dx 



Wa{x)dx 

= ^29+1 {XV) ■ / Wa{x)dx < OO. 

OO 

:= e-^^fc'^xVA;!, x > 0, 

k=l 

and V = Jj^al{x)dx is the volume of the ball -B2r, • D 

2.2.2. Percolation models. In Section [3] below, we explore the relationship 
between phase transitions in our model and two (related to each other) 
percolation models on the underlying graph G{'j,D). 

The continuum percolation consists in the appearance of an infinite con- 
nected component of G{'j,D), see |2H I22j . It is known that such a com- 
ponent exists with probability one (respectively, zero) if the intensity A of 
the Poisson measure Tr;^ satisfies inequality A > A* (respectively, A < A*), 
where A* > is a critical value, which depends on the connection function 
g. The threshold value satisfies inequality A* > with as in (j2.2p . 

Observe that there can only be a single infinite connected component, see 
Theorem 6.3, page 172]. 

The bond percolation model on G{'y,D) is formulated as follows. Let 
q £ (0, 1) be fixed. Then each edge of G(7, D) is marked independently open 
with probability q, and closed otherwise. Now we can form a new graph, by 
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removing closed edges, and consider the percolation problem thereon. To 
make this procedure consistent with the continuum percolation discussed 
above we introduce another system of random variables 

D := {Dj,y = Dyx, x,y eW^, x / y, \x - y\ < r*} 

on (ri, J^, P) such that each Dxy takes values 1 and with probability q and 
1 — respectively. All Dxy are mutually independent, are independent of the 
connection variables D and of the underlying Poisson point process. Let DD 
denote the system of product random variables {D^yDxy, x,y ^ M.'^, x ^ 
y, \x — y\ < r^,}. We say that G{'y,D) admits Bernoulli bond percolation if 
the graph G(7, has an infinite connected component. It is clear that 
the probability that given x,y £ 'y are connected in G{'~f,DD) is equal to 
qg{\x — y\). That is, the bond percolation in G(7, D) is equivalent to the 
continuum percolation with connection function qg. Hence, in accordance 
with the first part of this subsection, for a fixed intensity A > A* , there exists 

a critical probability G (0,1) such that the graph G{'j,DD) contains an 
infinite connected component with probability 1 (resp. 0) q > q.^ (resp. 
q < q^). It is clear that (7* > l/Xg*, cf. [7]. 

2.3. Tempered Gibbs states. In modern equilibrium statistical mechan- 
ics, the notion of thermodynamic phase of a system of bounded spins living 
on a fixed graph like Z'^ is attributed to extreme elements of the set of cor- 
responding Gibbs measures, see, e.g., [51 Chapter 7] or [Ml Chapter III]. 
However, for unbounded spins, not all extreme Gibbs measures may have 
physical meaning. It is believed that the measures corresponding to ob- 
served thermodynamic states should be supported on spin configurations 
with 'tempered growth' see e.g., [TH] or a more recent development in [31 dS] 
and [21 Chapter 3]. In this approach, only tempered Gibbs measures are 
taken into account, and hence a phase transition is related to the existence 
of multiple tempered Gibbs measures. Note that the theory of quantum sta- 
bilization and phase transitions in quantum anharmonic crystals developed 
in [H [3l [151 ES] and [21 Chapter 6] with the use of tempered Gibbs measures 
is consistent with the corresponding phenomena observed experimentally. 
Thus, in this work we take this approach and will study quenched Gibbs 
measures introduced in Definition 12.11 with a priori prescribed support prop- 
erties. We call them tempered Gibbs states. Clearly, we first have to show 
that such elements of Q('y,D) do exist, which is certainly not immediate 
if the underlying graph has unbounded vertex degrees and the single-spin 
distribution x has noncompact support. 
Thus, for an a > 0, we define 

(2.13) S(a) := |(j G M'^ : ^\ax\^Wa{x) < ool, 

where the weights Wa are as in (j2.9p . For each fixed 7, 5](a) is a Borel subset 
of and its elements are called tempered configurations. Then 

gt(7,^) := {/xGg(7,I)):MS(Q)) = l} 
is called the set of tempered (quenched) Gibbs states. 
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Theorem 2.3. Let the single-spin measure x be such that h2.3\) holds. Then 
the set of Gihhs states QtiljD) is almost surely nonempty. Moreover, for 
each positive 1!) and a, there exists an almost surely finite C{-d,a) > such 
that, uniformly for all fi S Gti'l, D), 

(2.14) / exp (^9 V|(J^|^^i'«(x) ) //(dcj) < C(^9,a). 

Proof. The statements fohow from Theorem 1 of |16j since all the condi- 
tions of that theorem are satisfied in view of ()2.6p . Proposition [221 and the 
assumed properties of x- D 

Let us make some comments. We first mention that the boundedness 
Jxy < (j)* assumed in (|2.5p and ()2.6p has been imposed for simplicity only. 
It can be relaxed by passing to 'tempered' interaction intensities as in [17j . 
Next, the existence of Gibbs measures follows from the relative weak com- 
pactness of the family {n^(-|cj) : A G Bh{^'^)}, for at least some a G 
which in turn follows from the tightness of this family by Prokhorov's theo- 
rem. A typical choice of a, for which the tightness is proven, is a^; = s G R 
for all a; G 7. Note that such a is tempered, see ()2.13p . Then the accumu- 
lation points of the family {n^(-|(j) : A G B\)(R'^)} are shown to obey the 
Dobrushin-Lanford-Ruelle equation and to satisfy the estimate in ()2.14p . 
in which the constant C(t?, a) can be expressed explicitly in terms of the 
weights as in ()2.9p and the parameters defined in (|2.10p . cf. Proposition 12.21 
From the latter fact one gets that the mentioned accumulation points are 
tempered, and hence belong to Gti'y, D). Note that Gti'y, D) is nonempty for 
all those u, for which both a{a, 9) and h{a) are finite. By similar arguments, 
one can show that Gtil-, D) is compact in the weak topology. 

Let {A„}„gN C i3b(M'^) be a cofinal sequence, which means that A„ C 
A„_l_i, n G N, and each A G B\:,{W^) is contained in a certain A„. The 
relative weak compactness of the family {n^(-|a) : A G i3b(IR'^)} yields 
that, for each s > 0, there exists a cofinal sequence {A„}„gN such that the 
sequence {11^ (•|o')}„gN with ax = s weakly converges to a certain element 
oigt{-i,D). For a > 0, by 

(2.15) ^±«gg,(^,D) 

we shall denote limiting elements of ^t(7,L)) which correspond to s = ±a. 
Note that each /x^*^ may depend on the sequence {A„,}„gN along which it 
was attained. 

Now we turn to the single-spin measure x- If it h^-s compact support, as 
was the case in [2^, then ()2.3p clearly holds for any u and x. The most 
known example of such x is 

(2.16) x{dt) = 6-i{dt) + 6+i{dt), 

which corresponds to an Ising magnet. Here 5s is the Dirac measure con- 
centrated at s G M. Since this magnet will be used as a reference model, 
we reserve a special notation ^^*^™§(7,Z)) for the set of all corresponding 
Gibbs measures, which are clearly tempered. By z/^ G ^^^^'"^(7,1?), we de- 
note the limiting Gibbs measures as in (j2.15p with o = 1. In this case, 
however, we have uniqueness, i.e., are independent of the sequences 
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{A„}„gN along which they were attained. This holds because, for each x 
andi/G^^^''^e(7,D), 

axV^{da) < / ax^idcr) < I axiy^{da). 

That is, ly'^ and are the maximum and minimum elements of ^^^^^^(7, D), 
respectively, cf. [HI Theorem 3.8]. 

In the case of 'unbounded' spins, a natural choice of the single-spin mea- 
sure is 

(2.17) x{dt) = exp{-V{t))dt, 

where y : R — )■ M is a measurable even function such that: (a) the set 
{t G M : V{t) < +00} is of positive Lebesgue measure; (b) V(t) increases at 
infinity as Itl""*"^ with some e > and u being as in (j2.3p . This includes the 
case where V is a polynomial of even degree at least 4 with positive leading 
coefficient, cf. [M El HSl HI] • 

2.4. The question of uniqueness. In the theory of Gibbs states, along 
with the existence issue the problem of uniqueness/nonuniqueness of such 
states is equaly important. Then a phase transition is understood as the 
possibility to pass from the uniqueness to nonuniqueness by changing the 
relevant parameters of the model. Thus, prior to proving almost sure non- 
uniqueness of n G Gtil , D) , which holds for large A and (p^, see Theorem 13. II 
below, we address the question of whether the same uniqueness does actually 
hold for some values of these parameters. For small enough A, see \22\ 
Theorem 1] or |21t Theorem 6.1, page 152], all the connected components of 
the graph G(7, D) are almost surely finite, and hence Gtil, D) is a singleton 
for all On the other hand, Qt{l, D) is a singleton if and only if, for each 
X € 7, arbitrary a G S(a), and any cofinal sequence {A„}ngN; one has 

(2.18) lim / ax^l {d(T\a) = 0. 

This equivalence holds for any symmetric ferromagnet satisfying the bound 
in ()2.14p . that can be proven by standard arguments based on the Strassen 
theorem, see [18] for more detail. Actually, for models with 'unbounded' 
spins living on an infinite connected graph, like the one we consider, there 
are no tool£| for proving (|2.18|) . For the Ising ferromagnet, the uniqueness 
in question on an infinite connected component of 6(7,!?) can be obtained 
by percolation arguments, see |9, Theorem 7.2]. For 'unbounded' spins, 
by means of the Brascamp-Lieb inequality for log-convex perturbations of 
Gaussian measures [5], we can prove the following result. If V in (|2.17|) is 
strictly convex and the energy functional has the form 

-El{aK\af^c) 

= ~2 Jxyi^fx ~ Cy) — — Jxy{<^x ~ O'y) ; 

{x,y}&E/^{-f,g) x&'y/^yedxn'YAc 

cf. (|2.4p . then, for all values of the intensity A of the underlying Poisson 
point process and all 0* > 0, the set Qt{t,D) is almost surely a singleton. 



-'^The celebrated Dobrushin uniqueness technique is not applicable here. 
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and hence no phase transition occurs. The existence in this case can be 
estabhshed exactly as in Theorem 12.31 above, cf. [6]. 



3.1. The statement. Recall that by a phase transition in the considered 
quenched ferromagnet we mean the fact that the set of tempered Gibbs 
states Qt{'j, D) almost surely contains at least two elements. It is equivalent 
to the appearance of a nonzero magnetization in states € Gt{l,D), cf. 
[HI Chapter 19] and (fZTHI) . 

Let us observe that there is no interaction between spins in different con- 
nected components of the underlying graph G{^,D). Then for a phase tran- 
sition to occur it is necessary that G{'y,D) almost surely possess an infinite 
connected component, that is, this graph admits a continuum percolation as 
described in subsection 12 . 2 . 2 1 ab ove . This is the case if the intensity A of the 
underlying Poisson point process obeys the bound A > A*. For A < A*, we 
have no infinite connected component in G(7, D) and thus \G{'y, D)\ = 1 with 
probability 1. In order to obtain a sufficient condition for a phase transition 
to occur, we will explore the well-known relationship between the Bernoulli 
bond percolation on the fixed sample graph G{'y,D) and the existence of 
multiple Gibbs states in the corresponding Ising model, established in [13], 
cf. subsection 12.2.21 above. Our goal is to prove the following result. 

Theorem 3.1. Let the measure x ^6 Theorem and such that 

x({0}) < x(I^)- Assume also that the intensity A of the underlying Poisson 
point process satisfies the condition A > A*. Then there exists a constant 
> such that, for any 4> satisfying ^2.6\) . the set Qt{'^,D) contains at 
least two elements with probability 1. 

The proof of this statement is based on the following result, cf. (|2.15|) . 

Lemma 3.2. Let the conditions of Theorem \ be satisfied. Then there 
exist a > 0, as in i2.15]) . and > such that, for cj) satisfying \2. b]) 
and for some o G j, the following estimate holds with probability one 



The proof of this lemma is given in the next subsection. 

Proof of Theorem \3.1[ - Since the integral in p.ip is the limit of those in ()2.18p 
with ax = a, then ()3.ip contradicts ()2.18p and hence implies non-uniqueness. 
On the other hand, by the invariance of x ^-nd of the interaction in (|2.4p 
with respect to the transformation a — )• — o" and cr — t- —a, we have 



Then (|3.ip yields /i"*"" ^ and hence the multiplicity in question. Note 
that o in (j3.ip belongs to the infinite connected component of G{'y,D), and 
the integral in ()3.ip is the mean value of the spin at this vertex in state /i"*"". 

3.2. Proof of Lemma 13.21 First, by means of percolation arguments of 
|13| . we prove the lemma in the case of the Ising model. Then we extend 
the proof to the general case by comparison inequalities. 



3. The Phase Transition 



(3.1) 
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3.2.1. The case of the Ising model. Recall that the single-spin measure of the 
Ising model is given in (|2.16p . ^^'^"§(7, D) denotes the set of all corresponding 
Gibbs measures, and iy+ G G^'"''''^ , D)) is the maximum Gibbs measure as 
in ()2.15p with a = 1. We are going to use the key fact proven in |13j : 
the Ising model with constant intensities Jxy = on the edges of an 

infinite Q, cf. ()2.4p . has at least two phases if and only if the graph admits 
the Bernoulli bond percolation with critical probability g* G (0, 1) such that 
0* > (log(l + g*) — log(l — q^:))/2. In our case, the graph G(7,D) with 
probability 1 admits this percolation and the threshold probability satisfies 
Q* > ^/^9*, cf. subsection 12.2.21 Then, for some o G 7, it follows that 

(3.2) / aou^ida) > 0, 

see Theorem 2.1] and also the proof of Lemma 4.2 therein. Here is 
the corresponding Gibbs measure of the Ising model with J^y = (/>*> 0. By 
the standard GKS inequality, see, e.g., \13\ Subsection 3.4], we have 

/ aoi'~^{da) > / aoi'~^{da), 
Jr-i Jr-i 

which together with p.2p yields the proof in this case. 



3.2.2. The general case. In this section, we estimate the integral in (j3.ip 
from below by the corresponding integral with respect to the maximum 
Gibbs measure 1^^ of the Ising model with a rescaled interaction intensity. 
Thus, from now on we explicitly indicate the dependence on <j) and use nota- 
tions Qtil, D, (j)) and ^^''^^(7, D, (f) for the corresponding sets of Gibbs mea- 
sures. The proof of the lemma immediately follows from the Wells inequality 
used, e.g., in [2l]. As the original publication [27] is hardly attainable, for 
the reader convenience we give a complete proof of this inequality here in 
the form suitable for our purposes. 

Proposition 3.3 (Wells inequality). Let a > be such that 

(3.3) x([aV2,+oo)) >x([0,a]). 

Then, for each x G 7 and each ^jl^"- G Q , D , cj)) , as well as for G 
gismg(^^2),a2,/.), we have that 

(3.4) / a^ix^^ida) >a [ a^v+ida). 

Proof. For the general choice of Ist IlJ^" be defined as in (|2.7p with 
Ox = +a for all x G 7. Each /x"*"" is the weak limit of {n^''^"}„gi!^ for some 
cofinal sequence {A„}„gN- In the case of unbounded spins, this convergence 
alone does not yet imply the convergence of the moments of n^'"*"" to that 
on the left-hand side of (j3.4p . Then we use the uniform in n bound as in 
(|2.14p . which can also be proven for all n^'"*"", and obtain 

/ a,nX'+"(d(T) ^ / a,;u+"(da), n ^ +00. 
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Since the sequence {A„}„gN is exhausting it contains a cofinal subsequence, 
{^nfelfceN, such that also 

/ cj^nj^''™^(dcr) / axV^{da), n +cx), 

where H^'^^"^^ is the kernel (|2.7p corresponding to the Ising single-spin mea- 
sure (|2.16p . interaction intensities a?J, and the choice Ux = +1 for all x G 7. 
Thus, the validity of (|3.4p will follow if we prove that, for each A which 
contains x, the following holds 

(3.5) / axTiT'id^) > « / <Txlf/''^Hda). 

Let Z/^{a) and Z^^^^{1) be the corresponding normalizing factors defined in 



. Then by (|2.7p we have, cf. 
(3.6) / a,.nX'+"(dcT)-a / axUl'''''>^{da) = (zA{a)z'f^^{r 

JRT JR-i ^ 

+ Y,{^y + aay]Ky\ (g) (xidax) ® x'^'"S(rf^^,)) ^ 

where x^'^^^^ is given in (|2.16p and Ky = a Ylz&dynA'= ■^y^ dy H A'^ 0, and 
Ky = otherwise. Note that Ky > in both cases. Then (|3.5p will follow 
from the positivity of the integral on the right-hand side of (|3.6p . Now we 
rewrite the integrand in (j3.6p in the variables := {ux ± adx )/V2, and 
then expand the exponent and write the integral as the sum of the products 
over X G 7A of 'one-site' integrals having the form 



{Sai / {utr'^{u-r^x{dax)0x''''''Hd^x) 
= Cx [ [{ax + ar^iax - a)"- + (a, - a)--(a, + a)"-] x{d(Tx), Cx > 0. 

JR 

Thus, to prove the statement we have to show that the integral on the right- 
hand side of (|3.7p is nonnegative for all values of mx,nx G Nq. By the 
assumed symmetry of this integral vanishes if rrix and Ux are of different 
parity. If both rux and Ux are even, then the positivity is immediate. Thus, 
it is left to consider the case where nix = 2k + 1 and Ux = 21 + 1. By the 
symmetry of X) it is enough to take k > I. Thus, we have to prove the 
positivity of the following integral 

/ \{a + af>^+\a - af+' + {a - af'+\a + af+'] x{da) 
Jr l J 



+00 



2/ (a^-aY^^ {a + af-^ + {a + a)'-' xido). 







k-l , /_ , „\k-l 



The function <^(o") := {a + a)'^ ^ + {a + a)'^ ' is increasing on [0, -|-oo). The 
integral on the right-hand side of the latter equality can be written in the 
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form 

r+oo 

(3.8) / - aY'^'^{a)x{da) = h{a) + h{a) + /3(a), 

JO 

Jo 

h{a) := r - a^f'+\{a)x{da) > 0, 

J a 

r+oo 

/3(a) := / (ct2 - a2)2'+V(a)x(d^) > a^'^'^{aV2)x{[aV2, +00)) 

In view of ()3.3p . the sum on the right-hand side of ()3.8p is nonnegative, 
which completes the proof. □ 
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